The algorithm for simulating generalized inverse gaussian variates is
copied from the R package fBasics from Diethelm Wuertz.
Dagpunar, J.S. (1989). An easily implemented generalised inverse
Gaussian generator. Communications in Statistics-Computation and Simulation,
18, 703--710.
Raible S. (2000). Levy Processes in Finance: Theory, Numerics
and Empirical Facts, PhD Thesis, University of Freiburg, Germany, 161 pages.